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For a dimcr with a non-degenerate orbital built from atomic wave functions of Gaussian shape we evaluate all 
the electron-phonon couplings derived from the one-body and two-body electronic interactions, considering both the 
adiabatic and extreme non-adiabatic limit. Not only the values of the coupling parameters in the two limits, but also 
the expressions of the corresponding terms in the Hamiltonian differ. Depending on the distance between the dimcr 
ions, some of the two-body couplings are comparable, or even larger than the one-body ones. 



I. THE MODEL. 



In a general two-site electron-phonon Hamiltonian H = H e i + H p h + H e i_ p h the interacting term H e i^ p h originates 
from developing H e j to first order in the ion displacements Ui (i = 1,2), where, in standard notation for a non- 
degenerate orbital: 
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We shall develop both the one-(e, t) and two-body(X, U, V, J z — J xy — P) electron interaction parameters, eval- 
uated as in Ref.l, by assuming a non-degenerate orbital described by Wannier functions built from atomic or- 
bitals of Gaussian shape. We associate to each site a Gaussian atomic-like orbital cj>i(r — Ri), with the ions 

centered at the positions R, = (±a/2 + u h 0, 0) (i = 1,2) . By defining N = (2/nf /4 T 3 / 2 , they read: 
4>i(r) = N cxp |— r 2 (x ± a/2 — Ui) 2 + y 2 + z 1 |. Then the Wannier functions 'J'i, ^2 can be written as: 

*i(r) = A(S)Mr) + B(S)Mr) * 2 (r) = S(5)0i(r) + A(5)0 2 (r) 



A(S) = [(1 + Sy 1 ' 2 + (1 - S)-^' 2 ]/2 B(S) = [(1 + S) 
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with S(u) = (0i |^>2) = exp[— r 2 (a + w) 2 /2], and u — u 2 —u\. To make clear our method of calculation, it is convenient 
to explicitate, as an example, the one-body local electronic energy: 



j ^*(r,R 1 ,R 2 ) 



2m 



V 2 + V 1 (r- fix) + V 2 (r-R 2 ) 



^ i (r,Ri,R 2 )d 3 r 



where the potentials originating from the ion cores at the displaced positions Ri and R 2 are: 

-1/2 



Vi = V(r - Ri) = -e 2 Z (x ± a/2 - Ul f + y 2 + z 2 



(i = l,2) 



(3) 



(4) 



with — e the electron charge, and -\-Ze the charge of the ion core. The local energy e (actually site-independent) can 
be decomposed into three terms, corresponding to the contributions from the kinetic operator (ey ) and from each 
one of the ionic potentials (e(/,ey , respectively). A similar decomposition holds for the hopping amplitude t. 
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We need to distinguish between the adiabatic and the non-adiabatic limit in evaluating the electron-phonon inter- 
actions, because the integrals have different kernels in the two cases. Indeed, in the adiabatic limit the displacements 
affect both the potentials V(r — Ri) and the Wannier functions \&j(r, R2) , expressing the requisite that the 
electronic charge distribution adjusts itself instantaneously at the position of the ions. We shall schematize the op- 
posite situation, where the electrons are slower than the ions, as realized by the electronic charge distribution staying 
centred around the undisplaced ion position, while the potentials are centred on the displaced ions. We shall call this 
the extreme anti-adiabatic limit. In the literature H the two limits are also named from, respectively, Frolich and 
Bloch. 



II. COUPLING TERMS IN THE ANTI-ADIABATIC LIMIT. 



In the anti-adiabatic limit e v does not change, therefore no electron-phonon coupling originates from it. The coupling 
terms derived from the two-body interactions are also identically vanishing in this limit, because they involve the 
Wannier functions and the inter-electronic Coulomb potential which are both insensitive to the displacements of the 
ions. The only non- vanishing electron-phonon non-adiabatic interactions arise from the variation of the potential 
contributions to e and t. Let us now succinctly describe their evaluation. Full details are given in Ref.3. The 

perturbation of e$ originates a term in the Hamiltonian connecting the local charge with the local deformation: 



S -(fo 1 ^ n io' u i + 9q Z >n 2aU2) ■ It has a formal similarity with the Holstein coupling term but it has a physically 
different origin, as Holstein [|J considered electrons moving along a chain of fixed spacing with vibrating diatomic 
molecules at its nodes. On site 1, its explicit evaluation H yields: 
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where F (x) — x 1 / 2 Erf(x 1 / 2 ). As, under site permutation, a — > —a, gjp = —g^\ as might have been anticipated 
by considering that, for equal charges n\ = ri2 and displacement amplitudes, with ei2 = (R2 — i?i)/a, the energies 

g^niUi • ei2 and g^n,2ii2 • e±2 on both sites coincide. Now symmetry requires u\ — —1*2 (a constraint which does 



not hold for the original Holstein model Qj) from which g^ = —gfr 
then reads: g$ '^2 a { n i a u i ~ n 2a u 2) The "crystal-field" coupling term g^' expresses the change in the energy e v , 
. To establish the form of this term, let us consider site 1, with charge n\. Its energy, after a displacement U2 of 
the ion on site 2, changes by an amount = g^P n i u 2 • e i2- This can be considered as the quantity measured 
by an observer sitting on ion 1 and watching the ion 2 moved by U2 . The equivalent measurement done by an 
observer on ion 2 watching the ion 1 displaced by u 1; yields E^ = g^pn2\i\ • e 12 . Assuming ri\ = rt 2 , andui = — u 2 
one must have E^ = E^ 2 \ implying g^p — ~ gff. Therefore for the dimer as a whole one writes this term as 
•9c/ Eu(' l iff u 2 — n2 a ui). The explicit evaluation j^] for site 1 yields: 



follows. The contribution to the Hamiltonian 
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From Eq.|6| and the change of sign of a under site permutation, g 



(12) 



The Su-Schricffcr-Hcegcr 
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characterizing the SSH Hamiltonian Hssh 



c\ a c lr7 )(u2 —Ux), is due to the modulation of the hopping amplitude t. To preserve the invariance 



of Hssh under site permutation, 7^ 12 ) = — ^ 21 ~> (see e.g refs.6 and 7). Indeed, its explicit expression M is: 
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Under site permutation a — ► —a and A — > B so that ^ 12 ^ = —jt 21 *) as expected. 

In conclusion, in the non-adiabatic limit the complete electron-phonon Hamiltonian is given by: 



H ep = 9a ^2(nicrUi - n 2< jU2) + gfp ^(nicr^ - n 2(T ui) 
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III. COUPLING TERMS IN THE ADIABATIC LIMIT. 



In the explicit expression of the different electronic interactions in the adiabatically displaced state, u invariably enters 
in the combination a + u. So.to obtain the corresponding couplings, one can simply take the derivative with respect 
to a of the parameters in Eqlfl as evaluated in Ref.l. 

There is some confusion inrhe literature about the correct form of the electron-phonon Hamiltonian obtained, in 
the adiabatic limit, from the variation of the local energy e, therefore we shall devote some space to clarify this point. 
In this limit, m ^ enters both the charge distributions and the potentials. The " crystal field" interaction couples 
the charge on site i to the relative position of site j through the modification of both the kinetic and the potential 
contributions to . 

Let's place the origin of the x-coordinate onto one of the displaced ions, at Ri, say. One has then to take into 
account the relative displacement of the ions. In the adiabatic limit, therefore, the overall e-derived electron-phonon 

coupling term in the Hamiltonian is: g^ 2 ^ ni a (u2 — Ui) + <?e 21 ^ ri2 a (ui — 1*2)- As g^ 12 ^ = —g^ 21 ^ we can write 
the total adiabatic contribution from local energy terms to the electron-phonon Hamiltonian as Hl p = g\ (ji2a + 

ni rT ){u2-'Ui). A Hamiltonian of this form was used in the papers of Ref.8, while those of Ref.9 proposed hamiltonianns 
incompatible with our results. After including the SSH term, the complete one-body electron-phonon Hamiltonian 
in the adiabatic limit has therefore the form |10| : 



K = 9i 12) £>la + ™2 CT )(«2 - «l) + 7 (12) + C lc la )(u 2 - Ul ). 
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Coming now to the explicit expressions of the one-body coupling parameters, let us write for conveniency g e 
£ 2) + & 2) and 7 (12) = 1 { V 12) + l { y 2) ■ We obtain: 
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Notice that, as the partial derivatives are linear in a, they change sign under site permutation. 

Also the two-body electronic interactions U, V, J(= P),X of Eq. (1) give rise to electron-phonon couplings, all of 
the form: 



rjep 
Y 



(12) 
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where Y = U,V,X,J and F(c\ a , c J(T ) is the function of Fermi operators representing the two-body interaction whose 

amplitude is Y. From the results of Ref.l, their evaluation is trivial. Notice that, as U(a) = J (a) + e 2 r/i/7r, then 
dU/da — dJ/da = dP/da. We list below their explicit expressions: 



dX 
da 



(-aTS) 



(1 + 3S 2 ) 
(l-S 2 ) 3 



1 + 2S 2 + F (a 2 r 2 ) - 2(1 + S 2 )F 



a 2 r 2 



— e 



— e 



S/a 



(1-S 2 ) 2 



4a 2 r 2 S' 2 



S/a 



du _ 2 r 

da 



7T 



(1-5 2 ) 2 

4 a r 2 S 2 ~ 



2(1 + S 2 



F 



S 2 - F (a 2 T 2 ) 

-Vs 



(1 - S 2 ) 3 
5 2 /a 



2 - S 2 + 2S 4 + S 2 F (a 2 r 2 ) - 4S 2 F 



a 2 T 2 



(1 -S 2 ) 2 



2a 2 r 2 



1 - 4S" - F (a 2 r 2 ) + 4F 



a 2 r 2 



(13) 



S 2 - F (a z T 



2 r 2\ 



F \^) ,N 



(14) 



d{V- J z /2) 
da 



4aT 2 S 2 



(1-S 2 ) 3 



3 - S 2 - 8S 4 - (7 - 5S 2 )F {a 2 T 2 ) - 4(1 - 3S 2 )F 



a 2 r^ 



1/a 



(1-S 2 ) 2 



2a 2 r 2 S' 2 



-1 -AS 2 + F (a 2 r 2 ) +4,Fo 



a 2 r 2 



25 2 + ^S 4 - f 2 - ^S* 2 ) F (a 2 r 2 ) + 25 2 



(15) 



Fig.l and 2 present the values of the one-body coupling constant for, respectively, the non- adiabatic, and the 
adiabatic limit, evaluated by assuming for the shape-controlling parameter of the Wannier functions T the typical 0] 

value l.A — 
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FIG. 1. Non 
r = l.oA -1 . 



■adiabatic coupling constants So ,Sc/ iT' U ' ( m e VA *) versus the dimer length (in A), evaluated assuming 
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The most unexpected result concerns g c f. While usually neglected in the literature M on metallic systems, this 
coupling has been recognized as relevant to polar materials jllj]. We find indeed that, when T = l.CjA -1 , g c f is larger 
than <7o for any a, and it becomes the largest parameter for a > 2.2A. For small a, the SSH coupling is the largest. 
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FIG. 2. Adiabatic coupling constants 
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Fig. 2 for the adiabatic case shows that g^ 12 ^ is always larger than the SSH interaction ^ 12 \ and particularly for 
large a there is an order of magnitude difference between them. 
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FIG. 3. Adiabatic coupling constants from two-body interactions dX/da, dU /da, d(V — J/2)/da (in eVA 1 ) versus the dimer 
length (in A), for F = 1.0 A" 1 . 



The case of T — 2. A" 1 (discussed in ref.3) shows that, the more localized are the orbitals, the more relevant is the 
role of g c t in relation to the other admissible couplings. 

Fig. 3 shows the couplings derived from the two-body electronic interactions for the same parameters as Fig. 2. In 
general, their values are much smaller than those of g e and 7^ 12 ) , with the possible exception of d(V — J/2) / da. Indeed, 

(i) 

that coupling arises from a physical mechanism not very different from the one originating gy , i.e. the vibration of 

the charge on site j as felt by site i. Similarly to g e also d(V — J z /2)/da decreases slowly with a, so that for large 
a those two are the only relevant couplings. Such interactions in the lattice have been recently discussed in ref.12, 
while their effects in the optical spectra have been treated in ref.13. 



IV. CONCLUSIONS 



We have presented the analytical evaluation of the electron-phonon coupling parameters derived from both one- 
and two-body electronic interactions in a model of a dimer with a non-degenerate orbital built from atomic orbital 
of Gaussian shape. We have shown that the coupling terms in the adiabatic and the anti-adiabatic limits differ 
qualitatively. 
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The evaluation of the coupling terms originating from the two-body electronic interactions shows that at least the one 
generated by the Coulomb repulsion between the charges on different sites, is comparable to, or even larger than, the 
couplings derived from the one-body interactions. The quantitative results for the coupling parameters, even if agreeing 
in order of magnitude with some estimates from experimental data |14fl are obviously model-depending. However, their 
ratios should be more close to the reality. In particular, the obtainea values of the coupling terms, when compared to 
the values of the electronic interactions resulting from the same Wannier functions Q| suggests that, for dimer lengths 
comparable to the lattice parameters in high temperature superconductors and colossal magnetoresistance materials, 
at most dll/da and dX/da can be safely dropped, while neglecting any of the other electron-phonon interactions is a 
questionable approximation. 
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